Nuclear profile dependence of elliptic fiow 
from a parton cascade 



Bin Zhang 

Department of Chemistry and Physics, Arkansas State University, 
P.O. Box 419, State University, Arkansas 724 67-0419, USA 

Abstract 

The transverse profile dependence of elliptic flow is studied in a parton cascade 
model. We compare results from the binary scaling profile to results from the 
wounded nucleon scaling profile. The impact parameter dependence of elliptic flow 
is shown to depend sensitively on the transverse profile of initial particles, however, 
if elliptic flow is plotted as a function of the relative multiplicity, the nuclear profile 
dependence disappears. The insensitivity was found previously in a hydro dynamical 
calculation. Our calculations indicate that the insensitivity is also valid with addi- 
tional viscous corrections. In addition, the minimum bias differential elliptic flow is 
demonstrated to be insensitive to the nuclear profile of the system. 
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Recently, the Relativistic Heavy Ion Collider has produced large amount of 
exciting new data. These new data give us valuable insight into the hot and 
dense nuclear matter. One of the important observables is elliptic flow which 
reflects the transverse anisotropy of particle momentum distribution. Elliptic 
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flow has been studied by many theoretical models, including non-Abelian en- 
ergy loss models [1,2], saturation models [3,4,5], parton recombination models 
[6,7,8,9,10,11,12,13], hydrodynamical models [14,15,16], and parton cascade 
models [17,18,19]. In this paper, we will study the elliptic flow using a parton 
cascade model. We will first introduce the elliptic fiow, and the parton cascade 
model used for this study. Then, we will use the parton cascade model to study 
the elliptic flow produced from two different initial transverse distributions, 
one proportional to the number of binary collisions, one proportional to the 
number of wounded nucleons. We demonstrate that even though, elliptic flow 
as a function of the impact parameter is very sensitive to the initial transverse 
distribution of particles, elliptic flow as a function of the relative multiplicity 
is almost independent of the transverse distribution. Furthermore, we show 
that the minimum bias differential elliptic flow is also insensitive to the initial 
transverse particle distribution. 

Elliptic flow is the elliptic deformation in the particle transverse momentum 
distribution [20]. It is usually characterized by the second Fourier coefficient 
of the particle azimuthal distribution [21]. If we use for the azimuthal 
distribution, and choose the azimuthal angle of the reaction plane to be zero, 
then 

f{(p) = vo + 2vi cos((/)) + 2v2 cos(2(/)) H . (1) 

The coefficient, V2, is the elliptic flow observable. It is the average of cos(20) 
of produced particles. If the transverse components of the momenta are know, 
it can be calculated by: 




(2) 



2 



In the above formula, (• • •) denotes the average over particles. Since initial mo- 
mentum distribution is isotropic, the elliptic flow is generated by flnal state 
interactions. Final state interactions (or pressure gradient) turn(s) the spatial 
anisotropy into momentum space anisotropy. It has been shown that the el- 
liptic flow is very sensitive to the initial stage evolution, and can be used as a 
sensitive probe of early dynamics [22,23]. 

In the following, we are going to use Zhang's Parton Cascade (ZPC) [24] to 
study elliptic flow at relativistic energies. The initial conditions are set up 
similar to those in a recent study by Molnar's Parton Cascade (MPC) [18]. 
In the local rest frame, the initial momentum distribution is thermal, with a 
temperature of 700 MeV. Particles are uniformly distributed between a space 
time rapidity of -5 and -|-5. The particle formation proper time is 0.1 fm/c. 
There are totally 2100 gluons per central event. As the momentum transport 
is determined by the momentum opacity, the following results are also correct 
if the total number of particles increases and the transport cross section de- 
creases by the same factor. To efficiently simulate momentum transport, we 
use isotropic differential cross sections that preserve the reaction plane of a 
collision. The total parton-parton elastic cross section agg will be varied to 
have values of 40 mb, 20 mb, 10 mb to study the response of the system. 
These cross sections are effective cross sections as no radiative energy loss or 
parton recombinations are included in the calculations. 

We will study the elliptic flow produced from two different initial transverse 
spatial distributions. One is proportional to the number of binary collisions 
per unit area. In this case, the particle number as a function of the impact 
parameter is also proportional to the number of binary collisions. The other 
distribution is proportional to the number of wounded nucleons per unit area 
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and the particle number as a function of the impact parameter is also propor- 
tional to the number of wounded nucleons. In generating the above distribu- 
tions, the nucleon-nucleon inelastic cross section crjviv is taken to be 40 mb, 
and the three parameter Woods-Saxon distribution is used for the nucleons 
inside one nucleus. These two spatial distributions are related to hard and soft 
particle production mechanisms, respectively [16]. 

The above initial spatial distribution and initial momentum distribution fac- 
torize. As local densities are sampled for the evolution of the expanding parton 
system according to the Boltzmann equation, the factorization is not automat- 
ically conserved. We also note that a geometry with sharp cylindrical nuclei 
always leads to larger elliptic flow values compared to the binary collision 
scaling case [1,17,18,25,26]. 
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Fig. 1. Number of particles as a function of the impact parameter for the binary 
collision (BC) scaling case and the wounded nucleon (WN) scaling case. Circles are 
generated by the simulation code. 

Fig. 1 gives the number of particles as a function of the impact parameter. 
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While the two distributions have the same number of particles when 6 = 0, the 
wounded nucleon scaling has more particles than the binary collision scaling 
case. In particular, when b = 10 fm, the wounded nucleon scaling has about 
twice as many particles as the binary coUision scaling case. Fig. 2 has the 
initial spatial ellipticity as a function of the impact parameter. The initial 
spatial ellipticity in the figure is defined through 

Note that this definition calculates the ratio first and then the average and 

the magnitude is smaller than the ratio of the averages. 
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Fig. 2. Initial spatial ellipticity as a function of the impact parameter. 

We first study the impact parameter dependence of the elliptic flow. The set 
up is similar to the recent MFC model study [18]. We approach the Boltz- 
mann limit by increasing the number of particles and at the same time de- 
creasing the cross section by the same factor [17,18,27,28]. In the binary col- 
lision scaling case, the rescaling factors for b = 0,2,4,6,8,10,12,14 fm are 
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A = 100, 100, 100, 220, 450, 1100, 5000, 50000. In the wounded nucleon scaling 
case, the rescahng factors are A = 100, 100, 100, 200, 300, 600, 2000, 11000. The 
convergence is checked by comparing to calculations done with A/2. The V2 is 
calculated for particles with a rapidity range of \y\ < 2. 

Prom Fig. 3, we observe that as the total cross section increases, or more pre- 
cisely, as the transport cross section increases, the elliptic flow increases. The 
binary collision scaling case is larger than the wounded nucleon scaling case 
for small impact parameters and smaller than the wounded nucleon scaling 
case for large impact parameters. This follows the trend of the initial spatial 
ellipticity. However, the V2 curves peak at smaller impact parameters than the 
e curves. This indicates that both initial ellipticity and initial particle density 
play roles in determining the elliptic flow. As the impact parameter increases, 
the elliptic flow increases with the initial ellipticity, however, after a point, 
the particle density is not high enough to generate enough response and the 
elliptic flow can not catch up with the initial ellipticity. It starts decreasing 
with increasing impact parameter. 

An alternative way of characterizing centrality is to use the relative central 
rapidity density, which is the ratio of the central rapidity density to that in 
central coUisions with 6 = 0. If we plot the elliptic flow as a function of the 
relative central rapidity density as in Fig. 4, we see that for the same transport 
cross section, the two curves corresponding to the binary collision scaling and 
the wounded nucleon scaling overlap. In other words, the impact parameter 
dependence of elliptic flow is cancelled by the impact parameter dependence 
of the multiplicity. Similar observations have been made in a recent hydrody- 
namics study [16]. This indicates that if we use the relative central rapidity 
density as a measure of centrality, the elliptic flow is not sensitive to whether 
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Fig. 3. Elliptic flow as a function of the impact parameter. Filled symbols are for the 
binary collision scaling case and open symbols are for the wounded nucleon scaling 
case. The curves are used to guide the eyes. Going from above, the three sets of 
results are for Ugg = 40, 20, 10 mb, respectively. For the binary collision scaling with 
Ggg = 40 mb case, the statistical error bars are also drawn. They are about the 
same for other curves. The diamonds are results for the binary collision case with 
Ggg = 40 mb and parton number rescaling factor of A/2. They agree well with the 
case with A particle division. 

the initial distribution is binary collision scaling or wounded nucleon scaling. 
It reflects the particle transport cross section, or in the case of hydrodynamics, 
the equation of state. Because of the viscous corrections, the cascade calcula- 
tions have a bend over when the relative central rapidity density is small. In 
contrast, the hydrodynamic calculations have an almost straight line depen- 
dence and overshoot data when the relative rapidity density is small. We also 
note that the Ugg — 40 mb binary collision scaling case is consistent with set 
D of [18]. 
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Fig. 4. Elliptic flow as a function of the relative rapidity density. Meanings of 
symbols are the same as those in Fig. 3. 

Now we turn to the study of pt differential elliptic flow which can give further 
information about the evolution [29]. Hydrodynamic studies agree well with 
low Pt data. At pt > 2 GeV, the data are consistent with a constant behav- 
ior while hydrodynamic results keep on increasing. In dynamic models, only 
when viscous effects are taken into account, is it possible to describe the devi- 
ation from the ideal hydrodynamical behavior. A recent hydrodynamic study 
demonstrates that the minimum bias pt differential flow is not sensitive to the 
initial nuclear profile. We want to know whether it is also true when viscous 
effects are taken into account. Fig. 5 shows the impact parameter averaged 
Pt differential flow. The calculations with binary collision scaling agree well 
with those with wounded nucleon scaling. This is true not only for the low pt 
region, but also for the high pt region where viscous effects are important. 



Another way of averaging over events is to calculate the multiplicity weighted 
average of cosine of the azimuthal angle. This gives the minimum bias elliptic 
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Fig. 5. Impact parameter averaged differential elliptic flow as a function of the 
transverse momentum. Meanings of symbols are the same as those in Fig. 3. 

flow. Results from the ZPC model are shown in Fig. 6. The binary collision 
scaling and the wounded nucleon scaling agree well with each other. This 
further demonstrates that the minimum bias pt differential elliptic flow is 
insensitive to the initial nuclear proflle. A comparison of Fig. 5 and Fig. 6 
shows that the minimum bias results are about 10% higher than the impact 
parameter averages. Hence the impact parameter averaged differential elliptic 
flow can be considered as a reasonably good approximation of the minimum 
bias differential elliptic flow. As pointed out in [18], the minimum bias elliptic 
flow weights in more central events. The central events in the wounded nucleon 
scaling case have lower elliptic flow than those in the binary collision scaling 
case. This can lead to a relatively smaller minimum bias elliptic flow in the 
wounded nucleon scaling case. However, the decrease in the relative amplitude 
can not be determined with the current statistics. 

In summary, we demonstrate that the elliptic flow as a function of relative 



9 



15 



minimum bias 




a =40,20,10 mb 
fined: BC scaling 
open: WN scaling 











1 



2 3 
p, (GeV/c) 



4 



5 



Fig. 6. Minimum bias differential elliptic flow as a function of the transverse mo- 
mentum. Meanings of symbols are the same as those in Fig. 3. 

central rapidity density in the binary scaling case agrees well with that in the 
wounded nucleon case. In addition, the minimum bias pt differential elliptic 
flow in the binary scaling case also agrees well with that in the wounded 
nucleon case. This is true not only for the low pt region, but also for the 
high Pt region where viscous effects are important and ideal hydrodynamics 
deviates from experimental data. As hadronization will not change the high pt 
elliptic flow [18], and recent research indicates the possibihty of extracting the 
parton elliptic flow from hadron elliptic flow [8] , the elliptic flow as a function 
of relative rapidity density and the minimum bias pt differential elliptic flow 
are promising observables for the extraction of information about final state 
interactions [30]. 
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